Submitted to the Astrophysical Journal, 201 1 September 4; accepted 2012 January 4 

Preprint typeset using J^TgX style emulateapj v. 5/2/1 1 



THE ATACAMA COSMOLOGY TELESCOPE: ACT-CL J0102-4915 "EL GORDO," A MASSIVE MERGING CLUSTER 

AT REDSHIFT 0.87 

Felipe Menanteau^^ John R Hughes \ Cristobal Sifon^, Matt Hilton^ Jorge Gonzalez^, Leopoldo Infante^, 
L. Felipe Barrientos^, Andrew J. Baker^ John R. Bond"^, Sudeep Das^'^'^, Mark J. DEVLIN^ Joanna Dunkley'^, 
Amir Hajian^, Adam D. Hincks^, Arthur Kosowsky^^, Danica MARSDEN^ Tobias A. Marriage^\ Kavilan Moodley^^, 
Michael D. NIEMACK^^ Michael R. Nolta^, Lyman A. Page^, Erik D. REESE^ Neelima Sehgal^"^, Jon Silvers^, 
David N. Spergel^, Suzanne T. Staggs^, Edward Wollack^^ 

^Rutgers University, Department of Physics & Astronomy, 136 Frelinghuysen Rd, Piscataway, NJ 08854, USA 
^Departamento de Astronomia y Astrofisica, Facultad de Fisica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
^School of Physics and Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 
"^Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, ON, Canada M5S 3H8 
^Berkeley Center for Cosmological Physics, LBL and Department of Physics, University of California, Berkeley, CA 94720, USA 
^Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ, 08544, USA 
^Joseph Henry Laboratories of Physics, Jadwin Hall, Princeton University, Princeton, NJ, 08544, USA 
^University of Pennsylvania, Physics and Astronomy, 209 South 33rd Street, Philadelphia, PA 19104, USA 
^Department of Astrophysics, Oxford University, Oxford, OXl 3RH, UK 
^^University of Pittsburgh, Physics & Astronomy Department, 100 Allen Hall, 3941 O'Hara Street, Pittsburgh, PA 15260, USA 
^^Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland 21218-2686, USA 
^^University of KwaZulu-Natal, Astrophysics & Cosmology Research Unit, School of Mathematical Sciences, Durban, 4041, South Africa. 
i^NIST Quantum Devices Group, 325 Broadway Mailcode 817.03, Boulder, CO, USA 80305 
I'^KIPAC, Stanford University, Stanford, California, 94305, USA and 
i^Code 553/665, NASA/Goddard Space Flight Center, Greenbelt, Maryland, 20771, USA 
Submitted to the Astrophysical Journal, 2011 September 4; accepted 2012 January 4 

ABSTRACT 

We present a detailed analysis from new multi-wavelength observations of the exceptional galaxy cluster 
ACT-CL J0102-4915, likely the most massive, hottest, most X-ray luminous and brightest Sunyaev-Zel'dovich 
(SZ) effect cluster known at redshifts greater than 0.6. The Atacama Cosmology Telescope (ACT) collaboration 
discovered ACT-CL JO 102-49 15 as the most significant Sunyaev-Zeldovich (SZ) decrement in a sky survey 
area of 755 square degrees. Our VLT/F0RS2 spectra of 89 member galaxies yield a cluster redshift, z = 
0.870, and velocity dispersion, (jgai = 1321 ± 106 km s~^ Our Chandra observations reveal a hot and X-ray 
luminous system with an integrated temperature of7i = 14.5±1.0 keV and 0.5-2.0 keV band luminosity of 
Lx = (2.19 ±0.11) X lO^^/z^Qcrg s~^ We obtain several statistically consistent cluster mass estimates; using 
empirical mass scaling relations with velocity dispersion. X-ray Yx, and integrated SZ distortion, we estimate 
a cluster mass of M2ooa = (2.16 ± 0.32) x 10^^ HjqMq. We constrain the stellar content of the cluster to be less 
than 1% of the total mass, using Spitzer IRAC and optical imaging. The Chandra and VLT/F0RS2 optical data 
also reveal that ACT-CL JO 102-49 15 is undergoing a major merger between components with a mass ratio of 
approximately 2 to 1. The X-ray data show significant temperature variations from a low of 6.6 ± 0.7 keV at 
the merging low-entropy, high-metallicity, cool core to a high of 22 ± 6 keV. We also see a wake in the X-ray 
surface brightness and deprojected gas density caused by the passage of one cluster through the other. Archival 
radio data at 843 MHz reveal diffuse radio emission that, if associated with the cluster, indicates the presence 
of an intense double radio relic, hosted by the highest redshift cluster yet. ACT-CL JO 102-49 15 is possibly 
a high-redshift analog of the famous Bullet Cluster. Such a massive cluster at this redshift is rare, although 
consistent with the standard ACDM cosmology in the lower part of its allowed mass range. Massive, high- 
redshift mergers like ACT-CL JO 102-49 15 are unlikely to be reproduced in the current generation of numerical 
N-body cosmological simulations. 

Subject headings: cosmology: observations — galaxy clusters: general — galaxies: clusters: individual (ACT- 
CL JO 102-49 15) — cosmic background radiation 



1. INTRODUCTION 

There are currently only a few examples of merging cluster 
systems in which there are spatial offsets (on the order of 200- 
300 kpc) between the peaks of the total and baryonic matter 
distributions. Some of these, like 1E0657-56 (the original 
"bullet" cluster at z = 0.296, Markevitch et al. 2002 ), Abell 
2146 (at z ■■ 



0.234, Ru ssell et al.||2010 ), and possibly Abell 
2744 (at z = 0.308, ,Merten et al.|pnTT| 



contain in addition 



^ Based on observations made with ESO Telescopes at the Paranal Ob- 
servatory under programme ID 086.A-0425. 



a cold, dense "bullet" of low entropy gas that is clearly as- 
sociated with the merger event. The offsets are due to the 
differing physical processes that act on the gas, galaxies, and 
dark matter. The gas behaves as a fluid, experiencing shocks, 
viscosity, and ram pressure, while the galaxies and dark mat- 
ter (or SO we posit) are collisionless. These bullet systems 
have been u sed to of fer direct evidence for the existence of 
dark matter ( [Clowe e t al. 2004, 2006) and to set constraints 
on t he self-interactio n cross-sectio n of the dark matter parti- 
cle ([Markevitch et al.|2004, RandaU et al.|2008HBradac et al.| 
2008| ). Additionally, the large merger velocity of 1E0657-56 
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of around 3000 km s"^ ( |Mastropietro & Burkert|2008t[Marke- 
lyitch 2006 ), required to explain the morphology and temper- 
ature of the gas from Chandra observations, is much higher 
than expected in cosmological simulations. These merging 
clusters provide interesting laboratories for several important 
topics in astrophysics and cosmology. 

Merging cluster systems can display a wide range of mass 
ratios, impact parameters, types of merging systems, and 
times since closest approach, and clusters evolve with cos- 
mic time. So it is worthwhile to find new bullet cluster sys- 
tems, especially at high redshift, to exploit their different ini- 
tial and final conditions and learn more about dark matter and 
the assembly of clusters through mergers. Massive clusters 
(~10^^ Mq) become increasingly rare at high z while their in 



fall speeds are typically higher in the early Universe (see |Lee 
|& Komatsu 2010). 

In this paper we repor t on ACT-CL JOl 02-49 15, a recently- 
discovered system (M arriage et al.pOll^ which appears to 
be an excellent new example of a bullet cluster at redshift 
z = 0.87. It was found by the Atacama Cosmology Tele- 
scope (ACT, [Fowler et al. 2007 ) through its str ong Sunyaev- 
Zel'dovich signal (SZ, Sunyaev & Ze ldovich|1 972 ) , an d was 
con firmed through optical and X-ray data by Menanteau et| 
|aL|(|20T0a). We demonstrate here that ACT-CL J0102-4915 



is a rare and exceptional system at a time in cosmic history 
when the Universe was only half its current age. Addition- 
ally, it is a contender for being the most massive and X-ray 
luminous galaxy clu ster at z > 0.6, likely competitors be- 

: = 0.89 with 



ing CL J1226+3332 ( |Maughan, et al.|pOQ4l ) at z = 0. 
masfjof M200C = (1.38 ±0.20) x 10'^ H^^Mq ( |Jee & 



& Tyson 

|200gr and SPT-CL J2106 -5844 at z ="l.l4 With ^^200^ = 



(1.27 ±0.21) X lO^^hj^MQ |Foley et al.|2"oTl] ). Therefore, to 
reflect its exceptional mass and to recognize the Chilean con- 
tribution to ACT we dub the cluster "El Gordo," which means 
The Fat/Big One in Spanish. 

Throughout this paper we quote cluster masses as M2ooa (or 
^500c) which corresponds to the mass enclosed within a radius 
where the overdensity is 200 (500) times the average (crit- 
ical) matter density, and we assume a standard flat ACDM 
cosmology with ftm = 0.27 and = 0.73, and give relevant 
quantities in terms of the Hubble parameter Hq = lOhjo km s~^ 
Mpc~^ . For this cosmology the angular diameter distance to 
the cluster is 1617 /z^q Mpc and V corresponds to a physical 
scale of 410hjQ kpc at the cluster redshift. Magnitudes are 
reported in the AB system and uncertainties are reported for 
68% confidence intervals. 

2. OBSERVATIONS 

We have embarked on an ambitious multi-wavelength 
follow-up program aimed at characterizing the baryonic (i.e., 
gas, stellar) and dark matter mass components and the SZ- 
mass scaling relations of the most massive clusters discov- 
ered by ACT via the SZ effect. As part of this effort we 
have obtained total mass estimates based on galaxy dynam- 
ics (via measuring spectroscopic redshifts with 8-m class tele- 
scopes, ^^^^^3^ 20 12), cluster X-ray emission with Chan- 
dra and stellar masses from infrared IRAC/Spitzer imaging. 
In the following sections we describe the different datasets 
and preview the results when applied to this study of ACT- 
CL J0102-4915. Images of the cluster from these various 
data sources are shown in Fig.[T] 

^ Throughout this paper we distinguish between masses quoted with re- 
spect to the average or critical density using the "a" or "c" subscripts. 



2.1. SZ observations 

ACT is a six-meter off-axis telescope designed for 
arcminute- scale millimeter- wave observations and it is situ- 
ated at an elevation of 5190 m on Cerro Toco in the Atacama 
desert in northern Chile (see [Fowler et al.|2010t [Swetz et aL] 



2011 for a complete description). One of the goals of this 



project is to measure the evolution of structure in the Uni 
verse using mas sive gala xy clusters detected via the SZ effect 
(e.g., S ehgal et al.|2011| ). 

ACT observations began in late 2007 and to date it has 
surveyed two sky areas: one near declination -55 degrees 
(the southern strip) and the other on the celestial equator. 
The 2008 southern observations covered a region of 455 deg^ 
and recently ACT has analyzed an additiona l 300 deg^ along 
the eq uator (centered on the SDSS Stripe82, |Abazajian et al. 
[2009] ). 

ACT-CL J0102-4915 was reported by ACT as a particu- 
larly strong SZ detection at a frequency of 148 GHz among 
23 high- significance c lusters from the 2008 southern survey 



( [Marriage et al.||20lT| . In Figure |T] (lower right panel) we 
show the filtered ACT 148 GHz mten 



intensity ma ps for ACT- 
CL J0102-4915 using a matched filtered (e.g., jHaehnelt & 

' - ' — — 



[Tegmark|19'96l|Melin, Bartlett, & D elabrouille"2606). 

The South Pole Telescope (SPT, Carlstrom et al. 2009 ) has 
reported its most significant cl uster detections on a 2,500 
deg^ survey d Williamson et al. 201 1| ), which overlaps with 
the ACT 2008 southern survey. This sample also contains 
ACT-CL J0102-4915; it is SPT's most significant SZ clus- 
ter detection to date by nearly a factor of two, with a com- 
parable beam-averaged SZ decrement to the Bullet cluster 
1E0657-56. Among all clusters in the combined ACT and 
SPT sample, ACT-CL JO 102-49 15 is exceptional in having 
the highest SZ signal (ACT, |Marria ge et al. 201 Ij) or the hig h- 
est significance SZ detection (SPT, Willia mson et al.|201 1| ). 

2.2. Optical Imaging 

The initial optical observations of ACT-CL J0102-4915 
were carried out during December 9-12, 2009 using the SOI 
camera on the 4.1-m SOAR Telescope in Cerro Pachon (09B- 
0355, PI: Menanteau) using the SDSS griz filter set with ex- 
posure times of 540 s (6 x 90 s), 720 s (6 x 120 s), 2200 s 
(8 X 275 s) and 2200 s (8 x 275s) in g, r, / and z respectively. 
These exposure times were twice the nominal observations for 
the program as ACT-CL JO 102-49 15 was one of the targeted 
detections suspected to be a higher-redshift cluster. The run 
conditions were optimal and all four nights were photometric 
with seeing in the range 0'/5 to Of' 9. A complete description 
of the observation s and data analysis is provided in Menan-] 
[teau et al.| ( [2010a| ). 

As part of our 201 OB spectroscopic follow-up campaign 
of the southern cluster sample, we secured additional deep 
and wider Riz imaging (intended for slit mask design) and ex- 
tensive multi-object spectroscopy (MOS) with F0RS2 on the 
Very Large Telescope (VLT, 086.A-0425 PLInfante). F0RS2 
is a multi-mode (imaging, polarimetry, long slit and multi- 
object spectroscopy) optical instrument mounted on the UTl 
(Antu) telescope. We devoted a total of 12 hours of observ- 
ing time (service mode) on VLT to ACT-CL J0102-4915 of 
which 10 hours were devoted to MOS and 2 hours to imaging. 

In our original confirmation of ACT-CL JO 102-49 15 we re- 
ported the larger concentration of galaxies to the NW as the 
optical cluster counterpart to the ACT source and used the 
brightest galaxy it contained as the cluster center (see Fig- 
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Figure 1. The multi- wavelength dataset for ACT-CL JO 102-49 15 with all panels showing the same sky region. (Upper left) The composite optical color image 
from the combined griz (SOAR/SOI) and Riz (VLT/FORS2) imaging with the overplotted Chandra X-ray surface brightness contours shown in white. The black 
and white inset image shows a remarkably strong lensing arc. (Upper right) The composite color image from the combination of the optical imaging from VLT 
and SOAR and IR from the Spitzer/IRAC 3.6/im and 4.5/im imaging. The overplotted linearly-spaced contours in white correspond to the matched-filtered ACT 
148 GHz intensity maps. Bottom left False color image of the Chandra X-ray emission with the same set of eleven log-spaced contours between 2.71 cts/arcsec^ 
and 0.03 cts/arcsec^ as in the panel above. The insert here shows the X-ray surface brightness in a cut across the "wake" region from the box region shown. 
Bottom right ACT 148 GHz intensity map with angular resolution of 1 M and match-filtered with a nominal galaxy cluster profile, in units of effective temperature 
difference from the mean. The color scale ranges from -85 /xK at the edges to -385 /iK at the center of the SZ minimum. In all panels the horizontal bar shows 
the scale of the image, where north is up and east is left. 

The F0RS2/VLT imaging aimed to obtain wider and deeper 
observations to confirm the substructure hint in the SOAR 
data and to secure galaxy positions for the spectroscopic 
observations. Our observations were carried out with the 
standard resolution of Of'25/pixel providing a field size of 
6/8 X 6/8 with exposure times of 640 s (8 x 80 s), 2200 s 
(10 X 220 s) and 2300 s (20 x 1 15 s) in R, I and z, respectively. 



ure[T]top panels). However, even in the smaller ^ 4/5 x 4/5 



coverage of the SOAR images (see Figure 9 from |Menanteau 
et al. 1 2010a), we noted a number of red galaxies with simi 



lar photometric redshifts trailing towards the SE. The lack of 
coverage and large photometric redshift uncertainties made it 
difficult to confidently identify these galaxies as part of the 
cluster. 
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Figure 2. Histogram showing the redshift distribution for the 89 cluster 
member galaxies and the Gaussian with the same mean and width (jga\. 

The optical data were processed using a slightly modified ver- 
sion of th e Python-based Rutgers Souther n Cosmology image 
pipeline ( jMenanteau et al.| [20091 POlObb that takes us from 
raw data to galaxy catalogs. Here we give a brief overview 
of the steps involved and refer the reader to the above pa- 
pers for a complete description. The initial standard CCD 
processing steps (i.e. flat fielding, bias subtraction) were per- 
formed via system calls to the ESO Recipe Execution tool 
(EsoRex]^ In order to increase the depth and coverage for 
ACT-CL JO 102-49 15 we co-added the imaging in the r, / and 
z and R, I and z from SOAR and F0RS2 respectively into 
an iiber dataset. As the F0RS2 RIz images were taken un- 
der clear but non-photometric conditions we utilize the SOAR 
photometric data to calibrate the photometry of the combined 
data. All science frames were aligned and median combined 
into a common pixel reference frame using SWarp (Berlin 
[][2006) to a plate scale of 0.''2525/pixel. Source detection 
and photometr y for the science catalogs were performed us- 
ing SExtractor ( [Berlin & Arnouts|1996| ) in dual-image mode 
in which sources were identified on the iiber /-band images 
using a 1.5(7 detection threshold, while magnitudes were ex- 
tracted at matching locations from all other bands. 

2.3. Optical Spectroscopy and Redshifts 

Based on the possibility that ACT-CL JO 102-49 15 was ex- 
ceptional, we devoted the bulk (10 of 15 hrs) of the awarded 
time on VLT in semester 201 OB to securing spectra for galax- 
ies with colors consistent with being early-type c luster mem- 
bers. We refer the reader to a companion paper ( jSifon et aL 
2012| ) where we fully describe the spectroscopic observing 
campaign of the 19 clusters at z > 0.3 in the ACT southern 
sample (of which ACT-CL J0102-4915 was part), including 
the data reduction, analysis and redshift determination. Here 
we summarize the most salient aspects of this work. 

The observations were executed in January 201 1 under pho- 
tometric conditions with seeing of < Of' 8. The observations 
cover the wavelength range ~ 4000 - 8000A using the GRIS 
300+11 Grism at 1.62 A/pixel and resolution A/ A A 660 at 
the central wavelength. At the cluster redshift ofz = 0.87, the 



[O II] 3727 A emission lines, the Call (K-H, 3950 A) (which 
is the spectral signature of elliptical galaxies), plus other ab- 
sorption lines such as the G-band (4300A) as well as the 
4000A break were located in the spectral range of the grism. 
These features facilitated the determination of galaxy red- 
shifts, galaxy cluster membership, and the galaxy velocity 
dispersion of the cluster. Integration times were 40 minutes 
per mask in order to maximize the number of galaxies while 
obtaining the necessary signal-to-noise ratio in the relevant 
lines. We observed a total of five F0RS2/MXU masks in this 
configuration, providing redshifts for a total of 123 objects of 
which 89 were cluster members. 

The redshifts were measured by cross-correlating the spec- 
tra with galaxy spectral templates of the Sloan Digital Sky 
Survey (SDSS) Data Release 7 ( Abazajian et al. 2009 ) using 
the RVSAO/XCSAO package for IRAF ( Kurtz & Mink 1998 ). 
Most of these member galaxies are ellipticals with no emis- 
sion lines, and only a few emission line galaxies have been 
found to belong to the cluster. 

Cluster membership was determined using a rest- frame cut 
in velocity space of 4000kms ~\ and then a pplying the shift- 
ing gapper method of Katge rt et al. (1996) in order to re- 
move galaxies that lie outside a given gap in velocity space 
from the main body of galaxies. In particular, we defined 
the main body in bins of 15 galaxies and used a velocity gap 
of 500km s~^ We used the well-accepted bi weight estimator 
( [Beers et ^[1990| ) to estimate both the mean redshift of the 
sample, z = 0.87008 ± 0.00010, and the velocity dispersior^of 
the cluster, cTgai = 1321 ± 106 km s~^. Errors were estimated 
via bootstrap resampling. The velocity dispersion of ACT- 
CL JO 102-49 15 is large r than aU other clusters in the ACT 
sample (Menan teau et al. 2010ai Sifon et al.[[2012| ) and also 
larger than SPT-CL 2106 -5844 which is the largest in the fuU 



2,500 deg2 SPT sample ( [Foley et al.|2"oTT 



^ http ://www.eso.org/sci/sof tware/cpl/esorex.html 



_ Williamson et al. 

2011). Using only the passive galaxies (i.e., objects with no 
emission lines) in ACT-CL J0102-4915, we obtained a mean 
redshift of z = 0.86963 ± 0.00012 with a velocity dispersion 
of 1280 ±108 km s-^ 

In Figure |2] we show the redshift distribution for all galaxy 
members from the F0RS2 spectra. We see no indication of 
substructure along the line of sight, confi rming the conclusion 
of Sifon et al. (2012) which used the Dressier & Shectman[ 
([1988j) test. However, as we will discuss in Sec. [3. 2[ we find 
substantial evidence for spatial substructure. 

2.4. Chandra X-ray Observations 

We were awarded 60 ks of Chandra observing time on this 
target during cycle 12 (PI: Hughes, Obsid: 12258) as part of 
our follow-up eff ort on the SZ clusters f ound in the ACT 2008 
455-deg^ survey ("Menanteau et al."2010a). The Chandra ob- 
servation was carried out on January 26-27, 2011 for an ef- 
fective exposure of 59,260 s using the ACIS-I array. Figure [T] 
(bottom left panel) shows the surface brightness of the cluster 
after point source removal, exposure correction, and adaptive 
kernel smoothing. 

The cometary appearance of the X-rays is remarkable, ex- 
tending even to the apparent presence of two "tails" extend- 
ing off toward the NW. Orthogonal slices of the X-ray surface 
brightness show that these "tails" arise from significant de- 
pressions in the X-ray intensity, at a level of 20-40%, in a 

^ As customary the quoted uncertainties for z and cTgai only reflect statisti- 
cal errors from 89 objects and do not include systematic effects. 
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Figure 3. Chandra X-ray spectra of ACT-CL JO 102-49 15 for the total clus- 
ter in black, the cool bright peak (region 1 in Fig.[5j in blue, and the highest 
temperature region (region 5) in red. The bottom panel shows the residuals 
between the best fit models (histograms in top panel) and data. 



100 



" 10 



+ 




Figure 4. The X-ray bolometric l uminosity vs . temp erature for a sample 
of well-studied clusters taken from Markevitch fl998V The Bullet cluster 
(1E0657-56) is the open square point at high temperature and luminosity 
from Markevitch. ^2006 > . ACT-CL J0102-4915 is the red circle in the same 
region of the plot. 

band extending from about V NW of the compact peak emis- 
sion region and continuing off to the NW out to the faint parts 
of the cluster. The depression is about 35^' (210 hjQ kpc) wide. 
We refer to this feature as the "wake". 

The integrated Chandra spectrum is shown in Figure [3] 
(top black points). The extraction region for this spectrum 
excludes the central peak emission (from within the inner- 
most contour in Fig.[T]) and extends out to the outermost con- 
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Figure 5. {Upper panel) Chandra X-ray image showing raw detected counts 
and the spectral extraction regions. {Lower panel) Best-fitted X-ray tempera- 
tures plotted versus spectral extraction region (see key in figure to left). Error 
bars are 1-cr. Region 1, corresponding to the bright peak, is the coolest region. 



Fig. [Tj which is at a radius of around 2!2 (1.0 /z^q 
An absorbed phabs^mekal model yielded a best- 



tour in 
Mpc). An 

fit (source frame) temperature of^r = 14.5±1.0 keV, metal 
abundance 0.19 ±0.09 with respect to solar values, and a 
bolometric luminosity of Lboi = 13.6 x IO^^^/z^q erg s~^ Fig- 
ure |4] shows the Lboi-7i relation with ACT-CL J0102-4915 
added as the red point (note that a redshift co rrection fac 
tor of (1 was applied to Lboi following Vikhlinin et 
2002| ). The nearby point on this figure is the closest 



al. 



comparison cluster, 1E0657-56. Similarly, the X-ray lumi- 
nosity of ACT-CL J0102-4915 in the 0.5-2.0 keV band is 
Lx = 2.19 ±0.11 X 10^5 erg s-^ 

The modest count rate of ACT-CL J0102-4915 (i.e., 0.2 s"^ 
corresponding to a total number of detected X-ray photons of 
nearly 12,000) limits the number of regions into which the 
cluster can be subdivided for spatially resolved spectroscopy. 
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Figure 6. The X-ray Chandra emission overlaid with the ACT 148 GHz data contours from Fig. fllf/^/fJ and the predicted SZ signal at 148 GHz from the 
deprojected X-ray data (right). The black cross represents the SZ position reported by WilHamson et aT. ( 201 1 1 for ACT-CL J0102-4915. The contour levels are 
linearly spaced between the limits shown in brackets on each panel. The minimum central value of the ACT and predicted signals agree to within ~15%. The 
white dotted line shows the symmetry axis chosen for the deconvolution. 



We perform this measurement by dividing the total emission 
into roughly equal parts as shown in Figure |5] Spectra were 
extracted and fitted to the same absorbed thermal model as 
before; the best-fit temperatures are plotted vs. region num- 
ber in Fig. [5] The coolest (region 1, ^r = 6.6±0.7 keV) 
and the hottest (region 5, kT = 22^^ keV) spectra are plot- 
ted in Figure^ in blue and red, respectively. Note the strong 
Fe K line (oDserved at E = 3.56 keV) in the spectrum of 
the coolest region, which yielded a rather high metal abun- 
dance of 0.57 ± 0.20. Elsewhere the Fe K line is not detected 
strongly, pointing to abundance variations in the cluster gas. 

We see that the peak emission region in ACT- 
CL J0102-4915 is an X-ray-bright, relatively cool, "bullet" 
of low entropy gas like in 1E0657-56. The enhanced Fe 
abundance suggests it was the core of a cool-core cluster 
that has, so far, survived a high velocity merger with another 
system. The cluster's high overall gas temperature and spatial 
temperature variations (Fig. [Sj argue for significant amounts 
of shock heating of the gas. The steep fall-off in the X-ray 
surface brightness toward the SE, as well as the "wake" in 
the main cluster gas toward the NW, clearly indicate that the 
bullet is moving toward the SE. 

The two lower panels of Fig. [T] show that the SZ and X- 
ray peaks are displaced by about l72. Similar offsets between 
the peak SZ and X-ray emission have also been reported in 
the bullet-like cluster Abell 2146 at z = 0.23, for which Chan- 
dra and Arcminute Microkelvin Imager (AMI) observations 
at 13.9-18.2 GHz show a remarkably similar separation of 
the SZ and X-ray signal distributi ons (see AMI Consortium:^ 
Rodriguez-Gonzalvez et al.|[2011| Figures 13 and 15). Two 
factors contribute to this. First, X-ray emission is more sensi- 
tive to denser parts of the cluster, since X-ray surface bright- 
ness is proportional to J n^T^dU where Ue and Te are the elec- 
tron number density and temperature and a ^i:^ 0.1 (for band- 
limited X-ray emission at high temperature), while the SZ sig- 
nal is proportional to the line-of-sight gas pressure / UeTedl. 
Second, the combination of Chandra and ACT's widely dif- 



ferent angular resolutions (< \" vs. \'A) is particularly rele- 
vant for observations of merging clusters with complex gas 
density substructure. 

In order to investigate these differences we developed a new 
technique for deprojecting X-ray images that, unlike tradi- 
tional methods utilizing azimuthally- averaged profiles, takes 
account of the full observed two-dimensional structure. The 
process assumes rotational symmetry about a user-defined 
axis and independently deprojects each column of pixels per- 
pendicular to that axis. For ACT-CL JO 102-49 15 we chose 
the symmetry axis so that it runs through both the peak X- 
ray emission at the bullet and the surface brightness depres- 
sion toward the NW. The axis is centered on R.A.=01 :02:54.9, 
Dec. =-49: 15:52.5 and oriented at a position angle of 136° 
(see Figure |6]). Beyond a certain radius from the center 
(1.'65 in the case of ACT-CL J0102-4915), where the sur- 
face brightness is background dominated, we used an analytic 
/3-profile fit to the azimuthally-averaged data to deproject the 
cluster emission smoothly to large radii. Interior to this ra- 
dius, using the adaptively- smoothed, 0.5-2.0 keV band im- 
age, we deproject the X-ray surface brightness in strips per- 
pendicular to the symmetry axis. Deprojection of the cluster 
emission is done in the usual manner ("peeling the onion") 
using the best-fit integrated spectral parameters given above 
to relate the X-ray emission to the gas density of each depro- 
jected shell. We deproject the north and south parts of the 
cluster separately, but interpolate the density value across the 
cluster when "peeling" back any individual layer (to avoid a 
discontinuity at the symmetry axis). In general this process is 
not guaranteed to produce positive definite values (as physi- 
cally required), so we force the computed density to be zero 
in cases where the deprojection would result in a negative (or 
more precisely, imaginary) value. We quote electron density 
values using an electron to proton number ratio of He /rip = 1.2 
and for mass densities we assume a mean mass of lAntp. 

The deprojection results in an estimate of the three- 
dimensional distribution of the gas density in the cluster. 
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From this we compute the predicted SZ signal at 148 GHz. 
Since the SZ depends linearly on the gas temperature, as a first 
approximation, we use the fitted X-ray temperatures in the 6 
regions shown in Fig. [3] to generate a temperature-corrected 
SZ map of the cluster. This was done by computing the ratio 
of expected SZ signal ( Nozawa et al.|199 8 ) for the observed 
temperatures in each of the six regions to the average cluster 
temperature and multiplying the single temperature SZ map 
by the appropriate ratio for the region. The multiplicative ra- 
tios vary from a value of 0.5 (region 1) to 1.6 (region 5). A 
more refined calculation would take into account the depro- 
jection of the X-ray temperatures, however, for our purposes 
here (i.e., to study the centroid differences), the simpler ap- 
proximation is sufficent. 

This temperature-corre cted SZ map is then convolved with 
the ACT beam model ( |Hincks et all [2010), added to a 
source-free region of the ACT map in the vicinity ACT- 
CL JO 102-49 15 (to include similar CMB and intrumental 
noise), and then is filtered with the same procedure used to 
detect ACT-CL J0102-4915. In Figure|6]we show contours of 
the predicted SZ signal overlaid on the Chandra data (right), 
next to an identical plot with contours of the ACT data (left). 
The resulting centroid difference between the peak of the SZ 
signal and the simulated SZ map is 6", which is within the 
expected positional error circle. Furthermore the intensity of 
the predicted SZ signal agrees to within ~15% with the ACT 
data. To conclude, we note that the offset between the peak 
X-ray and SZ signals is predominantly a result of the differ- 
ent sizes of the beams between Chandra and ACT (at least a 
factor of 60). 

2.5. Spitzer/IRAC Imaging 

As part of a larger warm-phase Spitzer program (PID: 
70149, PLMenanteau) which aims to characterize the stellar 
mass content of clusters in the ACT sample (Hilton et al., in 
prep.), in August 2010 we obtained IRAC 3.6/im and 4.5/im 
observations of ACT-CL JOl 02-49 15. Our observations were 
designed to provide coverage out to the virial radius for clus- 
ters at z > 0.4, using a 2 x 2 grid of IRAC pointings centered 
on the cluster position. A total of 10 x 100 s frames were 
obtained in each channel at each grid position, using a large 
scale cycling dither pattern. The basic calibrated data (BCD) 
images were corrected for pulldown using the softw are of 
Ashby & Horaj^and then mosaiced using MOPEX (Makovo? 
& Khan„200 5) to give maps which are ^ 13' on a side with 
0!'6 pixel scale. The maps for each channel were then reg- 
istered to a common pixel coordinate system. By inserting 
synthetic point sources, we estimate the 80% completeness 
depths of the final maps to be ^ 22.6 mag (AB) in both chan- 
nels. 

Matched aperture photometry was performed on the IRAC 
maps using SExtractor in dual image mode, using the 3.6/im 
channel as the detection band. We measure fluxes through A" 
diameter circular apertures, which are corrected to estimates 
of total magnitude using aperture corrections (as measured 
by [Barmby et al ||2008 ) of -0.35 ± 0.04, -0.37 ± 0.04 mag 
in the 3.6/im, 4.5/im channels respectively. The photometric 
uncertainties were scaled upwards by factors of 2.8, 2.6 in the 
3.6/im, 4.5/im channels respectively, in order to account for 
noise correlation between pixels introduced in the production 
of the mosaics which is not taken into account in the SEx- 



tractor error estimates. These scaling factors were d etermined 



using the method outlined in [Barmby et al.| ( [200^ . Finally, 
the uncertainties in the aperture corrections were added to the 
photometric errors in quadrature. 

3. ANALYSIS AND RESULTS 

3.1. Mass estimation 

Here we estimate the mass for ACT-CL JO 102-49 15 using 
all the mass proxies available from the X-ray data and optical 
spectroscopy. We present values in Table |2] and discuss them 
in detail in the next sections. In all cases we quote uncertain- 
ties that include measurement errors as well as the uncertain- 
ties in scaling relations between observables and cluster mass, 
relation. For the dynamical and X-ray masses, scaling to the 
common mass estimator M2ooa carries additional uncertainty 
that we quantify below and include in the Table. 

3.1.1. Dynamical Mass 

Here we report in more detail the dynamic al mass esti- 
mates for ACT-CL JO 102-49 15 calculated by |Sifon et al. 
(2012). The velocity dispersion of the cluster is dgai = 1321 zb 
106 km s~^ . This matches the largest reported dispersion of all 
known clusters at z > 0.6 f rom RX J0152.7-13 57 at z = 0.836 
and cTgai = 1322!^^ km s"^ ( [Girardi et al.|2005| ), although note 
that the quoted velocity dispersion is tor a well-studied merg- 
ing system which is resolved into two components of con- 
sider ably lower ma ss. We use the M2ooc-c^dm scaling relation 
from |Evrard et al.| (2008) to convert the measured velocity 
dispersion into a dynamical mass estimate. 



M200C = 



101 



C'"DM 



OJhjo(z) V^DM,15 



l/a 



Mo, 



(1) 



where hjoiz) = hoV^UX+zP + ^, (Jdmas = 1082.9 ± 
4.0 km s~\ a = 0.3361 ±0.0026 and (Jdm is the velocity dis- 
persion of the dark matter halo. The latter is related to the 
observed galaxy velocity dispersion by the velocity bias pa- 
rameter, by = cTgai/^DM- The latest physically motivated sim- 
ulations (see |Evrard et al . 2008, and references therein) indi- 
cate that galaxies are essentially unbiased tracers of the dark 
matter potential, (by) = 1.00 ±0.05. Using a bias factor of 
by = l for the velocity dispersion for all galaxies, we obtain a 
dynamical mass of M200fl,dyn = 1.86;^q 49 ^ IO^^/z^qM©. Simi- 
larly, using only the passively evolving galaxies without emis- 
sion lines, w e obtain M2ooa dvn = 1 .72+0^7 x IO^^/z^qM©. We 
note that the |Evrard et al.| ( |2008 |) scaling uses M200C as the 
mass within a radius where the overdensity is 200 times the 
critical density o f the Univers e, so the observed mass values 
reported in Sifo n et al. ( 2012| ) have been scaled from critical 



to average density using M2ooa = ^- ^^-nnt ^ ^looc - This con 



^ See http://irsa.ipac.caltech.edu/data/SPITZER /docs/| 

|dataanalysistools/tools/ contributed/irac/lixpulldowii/ 



version factor was derived using a [ Navarro, Frenk, & White 
(1997) mass profile (hereafter NEW) and the concentration- 
mass relation, c(M,z), from simulations (D uffy et al. 2008) 
at z = 0.87 for the mass of the cluster. The reported uncer- 
tainties in the conversion factor reflect the aiogc = 0.15 scatter 
in the log-normal probability distribution of c(M,z). The ra- 
dius r2oo = 21 1 1 ± 1 89 /z7Q kpc was also calculated using M2ooa 
and assuming spherical symmetry, using the average density 
of the Universe at redshift z, p(z) = +zf?>Hl/ (87rG). 

3.1.2. X-ray Mass 

We estimated the cluster mass from several X-ray mass 
proxies using the scaling relations in Vikhlinin et al. ( 2009 ) 
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and Kravtsov, Vikhlinin, & Nagai (20061. For the gas temper- 
ature and luminosity w e essentially lollow the prescriptions in 
Vikhlinin et al.| ( [2QQ9] ) to apply those scaling laws to the X-ray 
data. We obtain masses of M^mc = 1 •39!^ i^ ^ ^^^^ ^il ^0 
Msooc = 1.40!g:44 X 10^5 Mo from the Msooc-T^ and Msooc- 
Lx relations, respectively. We refrain, however, from apply- 
ing \hQ\l% "merging cluster mass boost" as implemented by 



Vikhlinin et al.| ( [2009 ), since, given the exceptional nature of 
ACT-CL JOl 02-49 15, it is not obvious that such a correc- 
tion factor is appropriate. Due to the complex substructure of 
the cluster, for the Mgas -based relations (Yx and Mgas), we use 
the 2D deprojection to generate an azimuthally-average radial 
gas mass profile about the cluster center, which agrees well 
(within ~15%) with the azimuthally- averaged core-excised 
surface brightness profile. For the Yx mass proxy we obtain a 
mass of Msooc = 1.55!^:}^ x lO^^/z^iM©. 

Since all of these mass proxies are potentially sensitive to 
merger boosts on Tx and L^, we also explored using Mgas, 
which is intrinsically unaffected by merger boosts. How- 
ever, the X-ray-derived measurement of Mgas, which comes 
from the reconstructed density by deprojecting the surface 
brightness, is sensitive to clumping. This can in principle 
result in a bias, depending on whether clusters are indeed 
more highly clumped during mergers. In the case of ACT- 
CL JO 102-49 15, the excellent agreement between our pre- 
dicted SZ signal from the Chandra data and the actual mm- 
band ACT data suggests that gas clumping is not significant. 
We use the scaling law for M^qqc -Mg as at redshift z = 0.6 from 
Kravtsov, Vikhlinin, & Nagai|(|2006), which yields a value of 



M500C = 1-67^0 20 X 1O^/z7qM0 and implies a gas mass frac- 
tion /gas = 0.133. The agreement among the four X-ray mass 
proxies is excellent: all are consistent within their \-cr ranges 
with a single mass value. 

Our M2ooa X-ray mass estimates are given in Table|2]using a 
conversion factor of 1.85;^q 2o ^^^^ ^sooc to M2ooa . This factor 
is derived using the same procedure as in Section [T.l.l[ 

3.1.3. SZMass 

We estimate the SZ-derived mass of this cluster using the 
peak Compton signal for th e brightest 0.' 5 pixel (}^7cmb) of 
490 ± 60/iK as pubHshed in ISehgal et al7| ( [20TT] ). We com- 
bined this measurement with a fiducial scaling relation be- 
tween jTcMB and M2ooa derived from simulations and cali- 
brated by X-ray observations with best-fit parameters given 
in that same work. Using this scaling relation, we find the 
mass of this cluster to be M2oo,sz = 1.64;^q 42 x IO^^/z^qM©. 
This mass is consistent with the indep endently-derived mass 
reported by SPT ( Williamso n"et al.|20 11 ) of M2oo,sz = (1-78 ± 
0.64) X 10^^ HjqMq, which uses a similar scaling relation be- 
tween the significance, and M2oo«- This value also includes 
a 6% reduction to the published mass value to reflect the 
change in redshift from their photometric estimate (z = 0.78) 
to the true spectroscopic value (z = 0.87). 

3.1.4. Stellar Mass from SED Fitting 

In order to measure stellar population parameters and the 
stellar mass content of the cluster using galaxy memb ers, 
we combine the optical photometry described i n Sec tion |2.2| 
with the IRAC photometry described in Section [23] by cross 
matching the catalogs using a If' 2 matching radius. We adopt 
the SExtractor MAG_AUTO measurements as estimates of 
total galaxy magnitudes in the optical bands, and the 4^' diam- 
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Figure 7. Examples of the SED fitting for four galaxy cluster members in 
ACT-CL JOl 02-49 15. The soHd line represents the best-fit SED for the 
broad-band g, r, i,z, 3.6/im and 4.5/im photometry shown as open symbols. 



eter aperture corrected magnitudes as measurements of total 
flux in the IRAC bands. 

We fit the broadband Spectral Energy Distributi ons (SEP) 
usi ng a si milar methodol ogy to that de scribed in Shapl ey et| 
laLl ([20051 see also Hilto n et al][2QT0l ). We augmentedlEe 
sample of 89 spectroscopic members to a to tal of 41 1 selected 
using photometric redshifts fBemt ez|2000| ) re-calibrated with 
the 5 17 r edshifts from seven clusters in the range 0.4 < z < 0.9 
from Sifon et al. (2012 ) and with the same griz photometry as 
ACT-CL J0102-4915. The improved photometric redshifts 
had typical & = 0.04 in the redshift range 0.8 < z < 0.9. We 
select galaxies with photometric redshifts within Az±0.06 of 
the cluster mean redshift. 

We construct a grid of Bruzual & Chariot] p00 3) solar 
metallicity models with exponentially declining star forma- 
tion histories with 20 values of r in the range 0.1-20 Gyr, 
and 53 ages in the range 0.001-7.0 Gyr. We adopt a ^Chabrier] 
( |2003| ) initial mass funct ion (IMF), and mode l the effect of 
dust extinction using the jCalzetti et al.| ( |2000| ) law, allowing 
E(B-V) to vary in the range 0.0-0.5 in steps of 0.02. The 
SEDs are fitted by analytically calculating the normalization 
N for each model SED where dx^/dN = 0, adopting the model 
with the lowest value as the best fit. The stellar mass 
is then estimated from the value of N. We therefore fit for a 
total of four parameters (age, r, E(B-V), M*). Errors on the 
SED parameters were estimated using Monte-Carlo simula- 
tions. Figure [7] shows some examples of fitted SEDs. 

Given the small number of SED points compared to the 
number of free parameters, coupled with degeneracies be- 
tween s everal parame ters (age, r, E(B-V); see the discus- 
sion in Shapley et al.||200 5 ), we do not expect the values of 
these parameters to be very robust. However, stellar mass is 
well constrained by the IRAC photometry, which probes the 
rest-frame near-IR at the cluster redshift, and stellar mass es- 
timates have been shown to be fairly insensitive to variations 
in other model p arameters, such as the ass umed star forma- 
tion history (e.g., |Forster Schreiber et al.|2004^ Shapley et al.| 
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2Q05| ). We note that uncertainties in the IMF and the modeling 
of thermally pulsating AGB stars (e.g., |Maraston|2QQ5||Con 



roy et al.|2QQ9| ) is likely to lead to the stellar mass estimates 



only being accurate to within a factor of two; these system- 
atic uncertainties are not taken into account in quoted errors 
on . The total stellar mass for all cluster members within 
r2oo = 2nihj^ kpc is M*oo = 1 .3 1 ± 0.26 x lO^^M© . This sug- 
gests a ratio of stellar mass to total mass of < 1% within r2oo 
for ACT-CL J0102-4915. 

3.2. Cluster Structure 

In the ACDM paradigm, the dark matter distribution should 
be aligned with the distribution of galaxies if the dark mat- 
ter particles are indeed collisionless. The cometary and 
bullet-like X-ray structure of ACT-CL J0102-4915 motivates 
our search for a spatial offset between the galaxies (that 
should primarily trace the dark matter distribution) and the 
X-ray gas in a cluster subcomponent as seen in, fo r example, 
IE 06 57-558 (Clowe et al. 2004 ), M S 1054-0321 (|Jeltemaet| 
|all2001| ), and MACS J 0025.4-1222 (|Bradac et al.|2008| ). 

Our VLT campaign ( Sifon et al.|2Q12| ) found no clear evi- 
dence for velocity substructure or a bimodal velocity distribu- 
tion in ACT-CL JO 102-49 15. However, visual inspection of 
the optical images immediately suggested an elongated struc- 
ture in the distribution of the red galaxies along the NW-SE 
axis, coincident with the direction of the structure seen in the 
X-ray emission (see Figure [T] left panels). We investigate the 
spatial cluster structure in more detail by constructing four 
maps: a) the surface density of galaxies, b) the /-band lu- 
minosity density, c) the 3.6/im/IRAC rest-frame luminosity 
density and d) the SED-fitted stellar mass density for spectro- 
scopically and photo metric ally selected cluster members as 



described in Section |3.1.4[ These four maps use a 6'7pixel 
grid smoothed with a Gaussian kernel of width 20'^ Fig- 
ure [8] shows our deep optical multi-band imaging of ACT- 
CL JO 102-49 15 overplotted with yellow contours of the four 
density maps (one per figure panel) along with the same X-ray 
surface brightness contours from Figure [T] displayed in white. 

In all cases the galaxy distribution is obviously double 
peaked and a bright cD-type galaxy sits near each peak. For 
the galaxy density (Figure [8] top left) and stellar mass den- 
sity (lower right) maps, the peak of the X-ray emission lies 
between the peaks of the density. The luminosity densities 
for the /- and 3.6/im-bands (Figure [8] upper right and lower 
left panels, respectively) are also double-peaked distributions. 
The /-band luminosity density, which probes the ultraviolet 
at z = 0.87, is somewhat unusual because it is dominated by 
the light contribution of the ultra-luminous Brightest Clus- 
ter Galaxy (BCG), discussed further below. As a result, the 
SE peak has the higher amplitude in the /-band, and coin- 
cides with the peak of the X-ray emission and with the lo- 
cation of the BCG. In the 3.6/im band, which better traces 
the stellar mass, the NW peak is dominant, as with the galaxy 
number distribution. Note that the spatial distribution of the 
spectroscopically-confirmed cluster members is also double- 
peaked, similar to the contours in the top left panel of Fig- 
ure [H 

Between and nearly equidistant from the two peaks in all 
four distributions lies a 20''-long straight blue arc, presum- 
ably the strongly lensed image of a high-redshift background 
galaxy (see black and white insert in Figure [TJ. This feature 
appears in the trough between the two peaks of all four distri- 
butions. 

Regardless of whether we inspect the individual galaxies 



or their integrated light contribution, it is clear that there is a 
double peak structure in their spatial distribution with nearly 
parallel contours that also align with the direction of the lens- 
ing arc. We use this feature to define two components in 
the cluster (which we refer to as the NW and SE compo- 
nents) and investigate their masses and differences in veloci- 
ties of the components. Specifically, we divide cluster galax- 
ies into two groups along the line defined by the pair of RA 
and Dec positions (01:03:22.0, -49:12:32.9) and (01:02:35.1, 
-49:18:09.8). The velocity dispersions of the NW and SE 
components are 1290 ± 134 and 1089 ±200 km s~\ respec- 
tively. We comp ute the ir masses using the same method de 
scribed in section 



3.1.1 



and obtain 1.76!g j2x 10^5 



H^qMq and 



1.06;^Q59 X IO^^/z^qMq, respectively. Formally this implies a 
range of possible mass ratios from approximately 5:1 to 1:1 
between the cluster components; in subsequent analysis we 
will consider a 2:1 mass ratio. Table [T] gives computed values 
for the cluster components, as well as the cluster total. The 
sum of the masses of the two components determined sepa- 
rately, (2.8 ±0.9) X lO^^hj^MQ, is about l-a higher than the 
total cluster mass inferred from the velocity dispersion of all 
galaxies. We discuss this more in the next section. 

The SE component has a line-of-sight peculiar velocity [^of 
586 ± 96 km s~^ with respect to the main (NW) component. 
We also note that the cluster BCG (z = 0.87014 ± 0.00030), 
which is spatially located in the smaller SE mass component, 
has a peculiar velocity of 731 ± 66 km s~^ with respect to the 
main (NW) component. In Figure [9] we compare the distribu- 
tion of the velocities for the NW and SE components shown 
in red and blue respectively against the distribution for the full 
sample. 

We also use the cluster spatial information to obtain the stel- 
lar mass of the components wit h their respective r2oo value 
using the SED fits from section |3.1.4| Following the same 
component segregration as described above, we obtain stellar 
masses of (7.5± 1.4) x lO^^M© and (5.6± 1.3) x lO^^M© for 
the NW and SE components respectively. 

3.3. Combined Mass 

Table E] compiles our mass estimates for ACT- 
CL 10102^4915 from X-ray, SZ, and optical data sets. 
While central values range over nearly a factor of two, 
all estimates are statistically consistent. This can be seen 
by noting that the 2-a mass ranges for all six mass esti- 
mates overlap for the mass range 1.78 x IO^^/z^qM© and 
2.88 X IO^^/z^qM©; in other words, for any mass in this 
range, none of the six mass proxies are more than 2-a away. 
For a combined mass estimate on which to base further 
conclusions, we construct a curve for a fit to a single 
combined mass value using the several mass estimates with 
their appropriate plus and minus error bars. For the three 
independent mass estimators from velocity dispersion, Yx, 
and Sunyaev-Zeldovich decrement, we obtain a minimum 
of 2.7 for 2 degrees of freedom (statistically acceptable) at a 
mass value of M2ooa = (2.16±0.32) x lO^^/z^^M©. We do not 
combine all six mass proxies (which gives a somewhat higher 
mass value), since the four derived from the same X-ray data 
are all likely correlated to a degree that is difficult to evaluate. 

The above error range includes uncertainties from measure- 
ment error, mass scaling relations, and extrapolation to a com- 

^ This is estimated as Av = c(zi -Z2)/(l the pecuHar velocity of the 
SE component in the reference frame of the NW component. 
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Table 1 

Dynamical Mass for ACT-CL JO 102-49 15 components 



Group 


iVgal 


redshift 


C^gal 

(km s-i) 


^200 

(kpc h-^) 


M200a 

(lOi^/z-iMo) 


ACT-CL JO 102-49 15 (NW) 
ACT-CL J0102-4915 (SE) 
ACT-CL J0102-4915 (total) 


51 

36 
89 


0.86849 ±0.00020 
0.87175 ±0.00024 
0.87008 ±0.00010 


1290 ±134 
1089 ±200 
1321 ±106 


2062 ±233 
1743 ±336 
2111±189 


1 76+0-62 




Figure 8. The isopleths from the joint spectroscopically and photometrically selected galaxies in ACT-CL J0102-4915 using four methods. The isopleth contours 
using (upper left) the number density of galaxies per square arcminute, (upper right) the rest-frame optical /-band luminosity, (lower left) the rest-frame 3.6/im 
Spitzer/IRAC luminosity and (lower right) the SED-fitted stellar mass. The luminosity contours are both in units of 10^ Lq arcsec"^, and the stellar mass map is 
in units of 10^ Mq arcsec"^. The X-ray emission is shown as white contours using the same levels from FigurefTl All panels agree on the bimodal structure of 
the galaxy distribution; moreover, both the galaxy number and stellar mass density maps have peaks with statistically significant offsets with respect to the X-ray 
emission peak. 
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Figure 9. Redshift distributions for spectroscopic member galaxies in the 
NW (red) and SE (blue) clusters subcomponent. The underlying distribution 
for all of the 89 cluster member galaxies for ACT-CL JO 102-49 15 is shown 
in gray. 



Table 2 

Mass measurements for ACT-CL J0102-4915. 



Mass Proxy 


Measurement 


Scaling Law M200a 






(1O15/Z-1M0) 



CTgal 


1321 ±106 km s-1 


M200 


-a 




Tx 


14.5 ± 1.0 keV 


M500C 


-Tx 


co+0.84 
^•-^^-0.59 


Lx 


2.19±0.11 


M500C 


-Lx 


2 60+1-36 


YX 


3.20 ±0.24 


M500C 


-Yx 


OO+0.78 
^•^^-0.55 


Mgas 


2.2±0.1 X IO^^Mq 


M50OC - 


Mgas 


3 10+0-92 


yT CMB 


490±60/iK 


M200 


-y 


1 64+0-62 



Combined 



2.16 ±0.32 



Note. — M500C masses from Lx and Tx were scaled to M2ooa 
masses assuming an NFW density profile and the mass-concentration 
relation of Duffy et al.( ( |200 8) as described in the text. Lx is in units 
of lO^^/z^Qcrg s~^ and is in the 0.5-2.0 keV band. Yx is in units of 
10^^ Mq keV. The combined mass is the minimum x of the mass prox- 
ies in bold (cTgai, Yx, and jTcmb). 



mon M2ooa mass scale. Estimates of galaxy cluster masses 
also may be impacted by a number of potential systematic er- 
rors. Since cosmological conclusions are strongly sensitive to 
the estimated cluster mass, it is important to have a realistic 
appraisal of these errors. Here we describe several relevant 
sources of systematic error or biases in both measured mass 
estimates and in the theoretical comparison with cosmological 
models. 

Dynamical mass estimates are based on the assumption that 
the system is dynamically relaxed, which is almost surely not 
the case in ACT-CL JOl 02-49 15, based on the observed spa- 
tial galaxy distribution. If we are observing a merger with the 



component velocities primarily in the plane of the sky, then 
significant kinetic energy can be contained in motions trans- 
verse to the line of sight which are not observed in galaxy 
redshifts; the dynamically inferred mass would be lower than 
the actual system mass. This seems to be consistent with the 
measurements of the dynamical masses of the two subcom- 
ponents in ACT-CL J0102-4915 (explored in Section |T2|), 



which when summed, (2.8 ±0.9) x IO^^/z^qM©, is in better 
agreement with the X-ray-derived masses than the full galaxy 
sample. Therefore the dynamical mass of the full cluster 
quoted in Table |2] may be an underestim ate. 

A number of groups (e.g., |Ricker & Sarazin 2001 ; Ran dall] 
et al.||20Q2| ) have shown that cluster mergers can result m a 



transient boost in X-ray temperaure and luminosity during a 
merger, compared to the X-ray emission of the merged clus- 
ter once it relaxes. The effect that such boosts have on the 
X-ray mass proxies has been more r ecently inve stigated by 
Kravtsov, Vikhlinin, & NagaiJ ( [200 61), 'Pool e et all ([20 07 ) and 
Yang et al.| ( |2010| ), among others. Vikhlini neraL] ^009 ) ar- 
gue that masses estimated from Tx should be increased by a 
factor of 1.17 (which we do not apply) for merger systems, to 
account for the fraction of merger kinetic energy that remains 
unthermalized. Kravtsov, Vikhli nin, & Nagai ( 2006) intro- 
duced the Yx mass proxy which has been shown to be fairly 



independent of merger dynamic s (e.g., [Kravtsov, Vikhlinin, &| 
[Nagai 2006; [Poole et al.|20Q7] ). ACT-CL J0102-4915 shows 
strong evidence for being in the merger process, particularly 
the compact region of X-ray enhancement and the double- 
peaked sky density distribution of galaxies. So in our com- 
bined mass estimate we strongly prefer to use Yx, but we note 
that the masses for all the X-ray proxies are consistent with 
each other. In particular the agreement of our masses from 
the Yx and Mgas proxies suggest that our X-ray masses are not 
significantly boosted by the merger. 

The accuracy of the mass scaling law for the SZ proxy (here 
yTcMB) is yet to be determined and in fact is one goal of the 
ACT follow-up program. Still it is interesting to note that 
the inferred SZ mass for ACT-CL J0102-4915 is only 60% 
of the 7x -inferred mass, albeit consistent within the large er- 
rors. The excellent agreement we find between the predicted 
SZ signal from the Chandra and ACT data (as illustrated in 
Figure [6]) suggests that we have a realistic gas model for the 
cluster and thus points to the mass scaling law, particularly at 
high redshift, as the main source of possible discrepancy. A 
similar difference between the SZ and X-ray derived masses 
is also seen for the massive cluster SPT-CL J2 106-5 844 at 



z = 1.132 (Foley et al. 2011), while Andersson et al. (2010) 



find that SZ-derived masses are some 22% less than their X- 
ray masses. Given six mass proxies for ACT-CL JO 102-49 15 
which are all statistically consistent, it seems unlikely that our 
combined mass estimate is subject to significant systematic 
error. 

For cosmological model comparisons, the theoretical mass 
function for clusters this massive is not precisely determined. 
In principle, direct numerical simulations of dark matter cos- 
mology in large volumes and low mass resolution can de- 
termine the high end of the mass function in a straightfor- 
ward way. The largest simulations to date have comoving 
volumes of roughly 2000 Gpc^, on the order of our Hub- 
ble volume, which nonetheless contain perhaps only a sin- 
gle cluster with mass 10^^ Mq at z = 1. An analytical model 
of the high-mass mass function accurate at the 10% level 
has not yet been demonstrated, so the mass functions used 
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Figure 10. The optical color-magnitude relations for ACT-CL JO 102-49 15 
for spectroscopically confirmed cluster members. Absorption and emission 
line galaxies are represented by red circle and blue squares respectively while 
the BCG is denoted by a yellow circle. The one spectroscopically confirmed 
AGN is represented by a cyan star. The dashed lines represent the fitted rela- 
tions for the cluster RXJ0152-1357 at z = 0.83 from Blakeslee et al. (2006) 
for the corresponding HST/ACS filters redshifted to z = 0.87. The gray tri- 
angle represents the observered magnitude of a M* galaxy (i = 22.16) in the 
CMR. 

for comparison involve extrapolation from lower-mass halos 
where statistics from numerical simulations are better. Care 
also must be taken to ensure that initial conditions in simula- 
tions are sufficiently accurate; the commonly used Zeldovich 
Approximation for initial conditions may underestimate the 
number of high-mass haloes by 10% to 30% compared to 
more accurate second-order Lagrangian perturbation theory 



Berlind, Farrar, &| 



([Crocce, Pueblas, & Scoccimarro||2006t 
Zaw|2010| ). Finally, large- volume simulations give dark mat- 
ter halo masses, which then must be translated into total halo 
masses, with possible systematic departures from the com- 
monly assumed universal baryon to dark matter ratio in large 
clusters. 

3.4. Cluster Evolution 

3.4.1. Color-Magnitude relations 

Considerable effort has been devoted to understanding the 
forma tion and evo lution of el lipticals in clu sters up to z — 1.4 
(e.g., Lidmanet al.| 2008 ; Hil ton et al.|2009| ). A consistent pic- 
ture has emerged in which these galaxies are a well-behaved 
class of objects with structural and chemical properties obey- 
ing simple power-law scaling relations. The high degree of 
homogeneity observed in these relations and in their stellar 



Figure 11. The optical to infrared color-magnitude diagram for cluster mem- 
ber galaxies. Symbols are as in Fig.Qo] The dashed lines ag ain represent the 
fitted r elations for the cluster RXJO 15^1 357 at z = 0.83 from |Blakeslee etal.l 
j2006) , now transformed to the Spitzer/IRAC bands as described in the text. 

population suggests formation at hi gh redshift followed by 



passive evolution of the ir galaxies ([B ower, Lucey, & Ellis 
1992; Barrientos, Schade, & Lopez-Cruz 1996; ElHs et al. 
1997 ; Holden et al. 2005 ; Mei et al. 2009 ). In principle, when 
precise photometry is available measurements of the intrinsic 
scatter about the color-magnitude relation can be used to con- 
strain the ages of the constituent stellar populations of cluster 
galaxies. 

Here we study the optical and optical/IR color-magnitude 
relations (CMR) for ACT-CL JO 102-49 15 to assess the evolu- 
tionary stage of the cluster. We construct the color-magnitude 
diagrams for spectroscopically confir med cluster members us- 
ing the basic approach d escribed in IBlakeslee et al. ( 2006 1 
and Menan teau et ar] ( |200 6 ). For total magnitudes we use the 
SExtractor MAG_AUTO values (with a Kron factor of 2.5) 
and we measure colors within the galaxy effecti ve radii 



which were computed using GALFIT (versio n 3.0, [Peng et al. 
2010| ) to fit each galaxy to a |Sersic| ( |1968| ) model from the 
PSF-convolved /-band images. Bright nearby objects are fit 
simultaneously while faint ones are masked. Magnitudes in 
each passband are then measured within a circular aperture of 
radius Re while enforcing a minimum radius of 1.5 pixel. 

In Figures [T0| and [TT] we show the optical and optical-to- 
infrared color-magnitu3e diagrams for spectroscopically con- 
firmed cluster members. Absorption and emission line galax- 
ies are represented by red circles and blue squares respectively 
while the BCG is denoted by a larger yellow circle. In order to 



ACT-CL J0102-4915, "El Gordo" 



13 



Look Back time [Gyr] 

5 6 7 8 




0.6 0.8 

Redshift 

Figure 12. The o bserved /-band magnitude s for the BCGs in the SZ cluster 
sample reported in Menanteau et al. ( 2010a ) compared to the BCGs of several 
high redshift clusters. We higWight the BCG of ACT-CL J0102-4915 as 
a yellow star. For comparison we show the curves for the expected z-band 
magnitude as a function of redshift for elliptical galaxies over a wide range 
intrinsic luminosities (i.e. the curves represent the observed magnitude for 
galaxies of the same luminosity). 

ease the comparison between panels all colors are shown with 
respect to the total /-band magnitude. In order to compare 
with the expected colors of passively evolving system in fig- 
ure 10 we show for comparison the fitted optica l CMRs for the 
cluster RXJ0152-1357 at z = 0.83 from Blakesl ee et al.| ( |2006l ) 
which was observed with the corresponding HST/A CS filter 
set and redshifted to z = 0.87 using a solar metallicity 'Br uzual 
& Chariot] ( |2003| ) burst model with r = 1.0 Gyr formeTat 
Zf = 5. We also display for reference the typical magnitude 
for an M* (/ = 22. 16) galaxy at the cluster redshift in the CMR 
as a gray triangle. In Figure [TT] we show the fitted i-z vs. / 
relation from Fig. ITQ| transformed to the Spitzer/IRAC band 
using the same stellar population model. 

Despite some scatter, a "red sequence" of passively evolved 
galaxies is quite distinguishable and it is in good agreement 
with the observed relations for RXJ0152-1357. However, un- 



like Blak eslee et al.| p006) we see bluer early-type galaxies 
brighter that M"" possiEly due to recent star-formation trig- 
gered by the merger of the cluster. We note however, that our 
classification of galaxies is only based on the color and visual 
information in our seeing-limited ground-based imaging. 

From the figures is also clear that the BCG is extremely 
luminous ('^1.5 mag brighter that any other galaxy) and sig- 
nificantly bluer in its optical colors than expected from a pas- 
sively evolving elliptical galaxy. In the following section we 
describe the BCG properties in more detail. 

3.4.2. Luminosity of Brightest Central Galaxy 

In Figure[T2l we plot the observed /-band magnitude of the 
BCGs of ACT clusters in the southern sample, plus the BCG's 
for three X-ray-selected, mass ive, high redshift clusters: 
RDCS 1252.9-2927 (Blakeslee et al.|2003t|Rosati et al.]20Q4l 
IDemarco et al ||2007|) at z = 1.237 XM MU J2235. 3-2557 
dMullis et al.| 12005 [ iRosati et al.||20Q9l) at z = 1.393 and 
X MMXC S J2215.9-1738 ( [Stanford eFari [20061 [Hilton et[ 
|aL[[2009]) at z = 1.457. For clarity the BCG for ACT- 
2-4915 is shown as the yellow star, with a luminosity 



of 22L* in the /-band. This is clearly an extremely lumi- 
nous object and it is the brightest BCG in the ACT southern 
sample. We also plot the expected apparent /-band magni- 
tudes of elliptical galaxies of different luminosities as a func 
tion of redshift. We use L* in the /-band as defined by 
[tonet al.^(2003) at z = 0. 1 (M* = -20.82+5 log(0.7) ) and 



passive evoluti on according to a solar metallicity, [Bruzual & 
Chariot (2003') r = 1 .0 Gyr burst model formed at z/ = 5. The 
resulting magnitude-redshift relation is plotted for a range of 
luminosities (0.4L*, L*, 4L*, 8L* and 16L*) aimed at repre- 
senting the brightest observed BCGs. 

The VLT/F0RS2 spectrum of the BCG in ACT- 
CL JO 102-49 15 has the signature of an E+A+[0 II] galaxy. 
This galaxy is a possible high-redshift analog of NGC 1275, 
the blue and radio-l oud cD galaxy in the Perseus cluster at 
z = 0.018 (McNama ra, Q'Connell, & Sarazin 1996). Blue 
BCGs like this one usually live in the centers of cool core 



ually 

clusters (e.g. Santos et al. [12011), and in this case the BCG 
is located within the cool, low entropy, merging core in ACT- 
CL JO 102-49 15 revealed by the Chandra data. We also obtain 
the physical size of the BCG from the str uctural parameters 
estimated using GALFIT in section [3AT] The galaxy has an 
effective radius Re = 10.4 hj^kpc for a Sersic index ^ = 1.3. It 
is likely that this BCG is the result of a recent galaxy merger; 
we speculate that this extreme galaxy resulted from a merger 
of two BCGs during the cluster collision. 

4. DISCUSSION 

4.1. A New Bullet Cluster? 

Given the absence of a weak-lensing mass reconstruc- 
tion, we use the galaxy distribution as a proxy for the to- 
tal mass distribution. This shoul d be a robust approxima- 
tion since for both 1E0657-56 ([Randall et al.[[2008j and 
MACS J005.4-1222 ( [Bradac et al.[[2008[ ) the galaxy and to- 
tal matter distribution closely track each other. Figure [8] (up- 
per left panel) shows that the SE galaxy distribution precedes 
the peak gas location by 17!' (173/z7Qkpc) in the approxi- 
mate direction of motion of the merger. This suggests ACT- 
CL JO 102-49 15 as another "bullet" cluster undergoing a high- 
velocity merger event resulting in a significant spatial sep- 
aration in the main baryonic (gas) and total mass compo- 
nents. Moreover, ACT-CL J0102-4915 is at a significantly 
earlier cosmic epoch, a lookback time around 4 Gyr longer 
than 1E0657-56. However, this picture requires confirmation 
through weak lensing mass mapping to determine the spatial 
distribution of the dark matter. 

4.2. Rarity of ACT-CL J0102-4915 

ACT-CL JO 102-49 15 is a rare, massive high-redshift clus- 
ter. We use the convenient fi tting formulas provided in ap- 
pendix A of [Mortonson et al. (2011) to compute the con- 
fidence level at which the existence of a cluster like ACT- 
CL JO 102-49 15 would be expected in standard ACDM cos- 
mological models, where bound objects grow via gravita- 
tional instability from Gaussian initial density fluctuations. In 
Figure [13] we show M2ooa(^) exclusion curves for which a sin- 
gle cluster with mass M200fl above the corresponding curve 
would conflict with ACDM and quintessence at 95% confi- 
dence level, including both sample and cosmological param- 
eter variance. In other words, the exclusion curves represent 
the mass threshold as a function of redshift for which a cluster 
is less than 5% likely to be found in a survey region for 95% 
of the ACDM parameter variance. 
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Figure 13. The exclusion curves, M200aiz), from the< Mortonson et al.H2011| 
fitting formulas where a single cluster above the correspondmg curve would 
conflict with Flat ACDM and quintessence at 95% C.L for both sample and 
parameters variance. The black dashed and dotted lines represent the exclu- 
sion curves for the full sky and the ACT survey analyzed region of 755 deg^ 
(i.e. the mass limit for a cluster at a given redshift that it is less than 5% 
likely to be found in 95% of the ACDM parameter probability distribution). 
We also show the same exclusion limit for 2800 deg^ (ACT+SPT area) as the 
black solid curve. The open circle shows our combined mass estimate for 
ACT-CL J0102-4915, with 1-cr error bars. 



For purposes of addressing how rare ACT-CL JO 102-49 15 
is, in Fig. [13] we plot the exclusion curves for the full sky, the 
region analyzed for the ACT survey (755 deg^) where ACT- 
CL JO 102-49 15 was discovered, and the combined ACT+SPT 
survey region (2800 deg^ of which 455 deg^ were observed 
by both ACT and SPT) as this is the most massive cluster at 
z > 0.5 in the combined area and detected by both experi- 
ments. Our combined estimated l-a mass range for the clus- 
ter (in TableB ranges from 1.84 to 2.48 x 10^^ h^^^MQ. At the 
cluster redsmit of 0.87, the lower end of this mass range falls 
below the ACT+SPT region exclusion curve, while the upper 
end of this range falls well below the full-sky exclusion curve. 
In other words, the cluster is not an unlikely occurrence in the 
ACT+SPT sky region, provided that its actual mass is l-a or 
more below our nominal mass of 2.16 x 10^^ HjqMq. If the 
cluster mass is 1 to 2-a above this nominal mass, then it is 
a rare occurrence in the ACT+SPT sky region but not unex- 
pected in the entire sky. The cluster mass would have to be 
about 3-cr above the nominal mass for it to lie above the 95% 
exclusion curve for the entire sky. We conclude that, while 
this cluster is clearly rare, it is not massive or early enough to 
put significant pressure by itself on the standard cosmological 
model. 

Because the high-mass cluster mass function is steeply 
falling with mass, the signals we observe are more likely 
to be from a cluster with mass below our nominal mass of 
2.16 X IO^^HjqMq and a positive measurement error than 
from a cluster with mass above this level and a negative mea- 
surement error. In other words, the cluster mass is more likely 
to be near t he lower end of our nominal mass range than the 
upper end. |Lima & Hu| ( [2005| ) give a discussion of this Ed- 
dington bias and compute a correction for this effect, which in 
our case is of order 5%. This correction does not qualitatively 
change any conclusions presented here. 



| Lee & Komatsu ( 2010| ) investigated the rarity of 
0657-558 in cosmological simulations down to z < 0.5 
based on its large merger velocity and found the existence 
of such a merging cluster to be in tension with the expecta- 
tions of ACDM. They suggested that such high merger speeds 
could potentially be more common at high redshift, while, on 
the other hand, massive, ~ IO^ ^Mq clusters are rare at ear- 
lier times. Using simulations, |Forero-Romero et aL] (|2010 1 



studied the expected probability distribution of the displace 
ments between the dark matter and gas cores, as observed in 
IE 0657-558, and found it to be expected in l%-2% of the 
clusters with masses larger than 10^^ Mq. 

Here we briefly consider the rarity of ACT-CL JO 102-49 15 
as a massive cluster undergoing a major merger with a mass 
ratio of ~ 2 to 1 at z = 0.87, ignoring its merger velocity and 
spatial separation between gas and dark matter. We exam- 
ined the output of a large- volumne (3.072 /z"^ Gpc)^ ACDM 
cosmological N-body simulation with 29 billion particles per- 
formed by I. Iliev with the Cubep3m code at NIC Juelich 
(G. Yepes, private communication). Within this large vol- 
ume there is only one cluster of mass 1.9 x 10^^ HjqMq at 
z — 1, which implies a very low number density of around 
10~^^Mpc~^ for such clusters. This system has only one small 
potential substructural feature with a mass of around 10^^ M©, 
suggesting a recent minor merger with mass ratio of 100 to 
1 . The comoving volume of the entire sky between z = and 
z = 1 is 54.7 h~^Gpc^ and hence the examined simulation sam- 
ples roughly half of the total cosmic volume accessible by SZ 
cluster observations. We conclude that a merger event of the 
mass and mass ratio similar to the one we are witnessing in 
ACT-CL JO 102-49 15 is a rare event in the current set of large 
volume simulations. Significantly larger simulations, or im- 
proved analytic approximations to halo growth and evolution, 
will be required to assess directly the likelihood of massive 
high-redshift merger events such as ACT-CL J0102-4915. 

4.3. Merger Speed 

We have argued above that the depression in the X-ray sur- 
face brightness northwest of the merging bullet is a "wake" 
caused by the passage of one cluster through the other. In 
order to study this feature in more detail, we examine the 
map of the deprojected electron density in the cluster mid- 
plane (i.e., the plane of the sky) of ACT-CL J0102-4915 (see 
Figure [14]) . The merging bullet is indeed dense with a peak 
density of = 0.047 cm~^, while the outskirts fall to a value 
He ^2x 10~^ cm~^. Formally our simple deprojection proce- 
dure fails in the region of depressed surface brightness (pre- 
sumably because the cluster is not axially symmetric as as- 
sumed) and the density values here were pegged to zero; this 
area appears black in Fig. [14] 

Nonetheless, clearly the NW region of depressed X-ray sur- 
face brightness requires a region of extremely low density in 
the cluster interior. The region of depressed density is about 
0.8-1' (370hjQ-460h^Q kpc) wide. It begins approximately 
2.6' (1200/z7^ kpc) behind the leading edge of the bullet to- 
ward the northwest and can be traced into the core of the clus- 
ter to within about V (460^0 kpc) from the bullet. At locations 
even closer to the bullet, the deprojected density map shows a 
complex, disorganized set of peaks and valleys suggesting the 
presence of merger-driven turbulence behind the cold merging 
core. 

The merger speed can be estimated from the peculiar veloc- 
ity difference between the SE and NW galaxy concentrations 
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Figure 14. Map of the electron density in the mid-plane of ACT- 
CL JO 102-49 15 from a deprojection of the Chandra image. The contours are 
at density values of 0.002, 0.003, 0.004, 0.006, 0.008, 0.011, 0.016, 0.023, 
0.032, 0.045 cm-3. 

as Vmerg = Vpec/ sin(6>) = 586/ sin(6>) km The collision is 
likely taking place close to the plane of the sky, otherwise the 
clear morphological features of that merger, i.e., sharp leading 
edge, post bullet turbulence, and the wake (Figure [14]), would 
not be so prominent. Working from this scenario we assume 
low inclination angles of ^ = 15° and 30° that provide merger 
speeds of 2300 km s~^ and 1200 km s~^. An accurate esti- 
mate of the merger speed will require detailed N-body/hydro 
simulations. 

4.4. Extended Radio Emmission around 
ACT-CL J0102-491 5? 



The re is increasing observational evidence (see Ferrari et al. 
2008[ for a review) that the large amount of energy released 
(^ 10^"^ erg) in massive cluster mergers can accelerate parti- 
cles to relativistic speeds that in the presence of an intracluster 
magnetic field can produce synchrotron radiation. These dif- 
fuse, non-thermal emission features, called radio relics or ha- 
los, have no obvious connection with individual cluster galax- 
ies and have typical sizes of about 1 Mpc. We searched for 
such radio features around ACT-CL JO 102-49 15 using Syd- 
ney University Molonglo Sky Survey (SUMSS) archival data 
at 843 MHz. The radio mosaic at the cluster location has an 
angular resolution of 45''x 6{)" (major axis aligned N-S) and 
a positional unce rtainty for bright sources (> 20 mJy beam~^) 
of l''-2'' (Mauch et al. 2003 ). In Figure [l5l (right panel) we 
show the SUMSS archival map with the location of three ra- 
dio sources, labelled A, B and C, indicated with blue crosses. 
Only source A was cataloged in Mauch et al.| ( |2003 ); the lat- 
ter two sources (B and C) were analyzed by us using IMFIT 
within the Astronomical Image Processing System (AIPS). 
Under the assumption that the typical position error for faint 
sources scales as the beam size times the signal-to-noise ra- 
tio, we estimate positional errors for sources B and C of ~ 1" , 
which we use in the search for possible counterparts (galax- 
ies, AGN) in our optical and X-ray data. There is an excel- 
lent compact optical, IR, and X-ray counterpart for source B 



(< 1^' away), making an AGN identification for this source 
likely. The nearest X-ray source to position C is 8^' away 
which is at the limit of the radio source positional uncertainty 
so a clear identification is not possible. Furthermore this radio 
source is located right at the apex of the merging bullet, at the 
plausible location of a bow shock. There is no plausible X- 
ray counterpart for source A, which is also extended approx- 
imately in the direction perpendicular to the merger axis and 
centered near where the countershock should appear. More- 
over, Figure p?] (left panel), an unsharp-masked image of the 
Chandra data~(the difference of two images convolved with 
2D Gaussian smoothing widths of a = 7.5^' and 30^0. reveals 
high frequency structure, i.e., a sharp edge to the X-ray sur- 
face brightness, at the locations of both radio sources A and 
C, lending credence to the idea that we are seeing radio relics 
similar to t hose seen in s everal nearby clusters, su ch as Abe ll 
3667 (Rottgering et al. 1997 ), Abell 3376 (Bagchi et al. 2006 ), 
and Abell 1240 ( Kempner& Sarazin 2001) that are beheved 
to be powered b y merger shoc ks (see, e.g., [van Weeren et al. 
[201 Ij ). jKempner et al.| ( |2004| ) suggest the term radio gischt 
for these features to distinguish them from other types of ra- 
dio sources in clusters, such as AGN relics and radio phoenix 
sources. 

The monochromatic radio power of source A, assuming 
it is located at the redshift of the cluster, is P14 = 7 x 
10^^ W Hz~^ Likewise the radio power of source C is 
Pi4 = 3 X 10^^ W Hz~^ The rest- frame frequency of the 
843 MHz observations at the redshift of ACT-CL JO 102-49 15 
is 1.576 GHz, which we multiply by a factor of 1.126 to con- 
vert to the monochromatic power at rest-frame frequency 1 .4 
GHz, which assumes a spectral index of -1. The radio power 
of source A is one of the most intense known among radio gis- 
cht sources. Feretti ( 2002 ) show that the radio power of gischt 
relics correlate with the cluster bolometric X-ray luminosity 
(albeit with large scatter). The high radio power and bolomet- 
ric X-ray luminosity of ACT-CL JO 102-49 15 put it at the ex- 
treme upper end of this correlation. For ACT-CL JO 102-49 15, 
an important caveat of our interpretation of the SUMSS mea- 
surements is that the angular resolution of the radio images 
is insufficient to identify and remove radio source contamina- 
tion, so the values of radio power quoted above for sources A 
and C should be considered upper limits to the diffuse radio 
emission. Approved observations at the Australia Telescope 
Compact Array (ATCA) will provide higher resolution imag- 
ing that will help to remove much of this uncertainty due to 
unresolved radio source emission at the site of the relics. If 
confirmed this will be the most distant radio relic observed. 

5. SUMMARY AND CONCLUSIONS 

We present a multi-wavelength analysis of ACT- 
CL J0102-4915 "El Gordo," a new SZ cluster which 
appears to be the first reported high-redshift Bullet-like clus- 
ter undergoing a major merger. We have used a combination 
of optical (VLT), X-ray (Chandra), and infrared (Spitzer) 
data to investigate the cluster's physical properties and find 
that it is an exceptionally massive. X-ray luminous and hot 
system, resulting in the highest SZ signal in the ACT survey. 
From the optical spectra of 89 galaxies confirmed as cluster 
members, we obtain a spectroscopic redshift (z = 0.87008) 
and velocity dispersion (cTgai = 1321±106 km s~^) for this 
cluster. Our 60-ks Chandra observation on ACIS-I provides 
a spectroscopic temperature of Tx = 14.5 ± l.O keV and 
bolometric luminosity of Lboi = 13.6 x IO^^/z^q erg s~^ From 
mass scaling relations for the velocity dispersion. X-ray Yx, 
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Table 3 

Summary of SUMSS radio source positions and fluxes near ACT-CL J0102-4915. 



Source 


Name 


RA 


Dec. 




size 


RA. 


in SUMSS 


A 


SUMSS J010246-491435 


01:02:46.8 


-49:14:35.7 


18.2±1.5 


70''x 59^' 


64° 


yes 


B 


SUMSS J010252-491316 


01:02:52.4 


-49:13:16.7 


7.2±1.2 


unresolved 




no 


C 


SUMSS J010301-491705 


01:03:01.1 


-49:17:05.0 


7.8±1.2 


unresolved 




no 



Note. — Only the first object (A) was cataloged, the latter 2 sources were analyzed by us using IMFIT (AIPS) on 
the mosaic image. Sources are ordered in increasing right ascension. 




Figure 15. Unsharp -masked Chandra 0.5-2.0 keV image for ACT-CL JO 102-49 15 created from the subtraction of two images convolved with Gaussians of 
a = 1".5 and 30^^ with the high-frequency features highlighted as thick white curves (left). The SUMSS 843 MHz image of the same region showing the position 
of the radio sources A, B, and C on the map and their respective contours level are shown in red {right). In both images and for reference we overlay the same 
X-ray contours used in Figure[5] 



and SZ distortion, we determine the total mass of the cluster 
to be M2oo« = (2.16 ±0.32) x lO^^/z^^M©, estabHshing it as 
the most massive cluster known at z > 0.6. 

From our warm-phase Spitzer IRAC and deep optical imag- 
ing, we measure the total stellar mass and constrain the stel- 
lar content of the cluster to be <1% of the total mass, in 
broad agreement with other massive clusters. The optical {riz 
bands) and IR (3.6/im, 4.5/im) color-magnitude diagrams for 
spectroscopically-confirmed members show a clearly-defined 
red sequence consistent with passively evolving galaxies, in- 
cluding a small fraction of actively star-forming galaxies. 
However the BCG is extremely luminous, significantly bluer 
than the red sequence, and appears as an E+A+[0 II] galaxy 
from the VLT spectra. Such blue BCGs are normally seen 
in the centers of cool core clusters and indeed this BCG sits 
right within the cool, low entropy, enhanced-metal-abundance 
merging core in ACT-CL JO 102-49 15 revealed by the Chan- 
dra data. 

While clusters as massive as ACT-CL JO 102-49 15 are quite 
rare at its redshift, the cluster does not pose any signifi- 
cant difficulty for the standard ACDM cosmology provided 
its mass is in the lower portion of its statistically allowed 
mass range. Our Chandra and VLT observations addition- 
ally show that ACT-CL JO 102-49 15 is undergoing a ma- 



jor merger with a mass ratio of approximately 2 to 1 be- 
tween its subcomponents. We find no analogous high-mass 
merging systems, with properties broadly similar to ACT- 
CL JOl 02-49 15, within any current large- volume cosmolog- 
ical N-body simulations (e.g., MICE, CubepSm). Thus the 
high mass and merger ratio at this redshift argues for an ex- 
ceptional rarity of ACT-CL JO 102-49 15 within our Universe. 
We expect that more detailed analysis of large cosmological 
simulations and better understanding of the propagation of 
the initial Gaussian fluctuations, particularly at high redshift, 
will be required to compare with the predictions of a ACDM 
model using massive systems like ACT-CL JOl 02-49 15. 

The X-ray surface brightness of our Chandra observation 
reveals a "wake" in the hot gas distribution that we attribute to 
the recent passage of one cluster through the other. The depro- 
jected gas density distribution shows that the wake requires a 
low-density cavity in the interior of the cluster. Understand- 
ing the evolution of these features, determining the merger 
speed and quantifying the rarity of a merging cluster system 
with the relative masses and velocities of the subcomponents 
of ACT-CL JO 102-49 15 will require further simulations and 
analysis. 

Inspection of archival SUMSS radio observations at 
843 MHz reveals the presence of radio sources located at 
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the southeast and northwest edges of the cluster, which we 
tentatively identify as radio gischt relics hosted by ACT- 
CL J0102-4915. They lack obvious X-ray counterparts and 
their positions relative to the cluster are consistent with the 
expected locations of merger-induced shocks. The northwest- 
ern radio source is intense (Pi. 4 = 7 x 10^^ W Hz~^) and is 
extended perpendicular to the merger direction. If upcoming 
ATCA observations confirm the diffuse nature of the radio 
emission, then ACT-CL JO 102-49 15 will be the highest red- 
shift cluster with a radio gischt relic. 

The next most critical observational step is to obtain strong 
and weak lensing mass estimates and maps from the ground 
and with the Hubble Space Telescope, which we are pursu- 
ing. This will allow us to confirm the presence of an offset 
between the baryonic and dark matter distributions, as well as 
provide a robust mass estimate unbiased by merger state or 
cluster astrophysics. A deeper approved Chandra observation 
will allow better mapping of the spatial temperature distribu- 
tion in the cluster, which may allow identification of shocks 
in the gas. It is now possible to use ALMA to search for these 
same features directly in the SZ effect using the low frequency 
channel (Band 3 covering 84-116 GHz). Confirmation of the 
radio relic is clearly a high priority since these are typically 
found in merging clusters; moreover ACT-CL JO 102-491 5 's 
high redshift will allow investigation of their cosmic evolu- 
tion. As we have outlined above, several theoretical steps also 
need to be pursued in order to fully understand the many dis- 
tinctive features of ACT-CL JO 102-49 15 that we have pre- 
sented here. We hope that this new exceptional cluster will 
fuel more in-depth simulations of cluster mergers and moti- 
vate continued searches for similar high-redshift systems. 
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